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SCATTER AND BEAM HARDENING 
CORRECTION IN COMPUTED TOMOGRAPHY 
APPLICATIONS 

BACKGROUND 

[0001] The invention relates generally to industrial computed tomography (CT) 
systems, and more specifically to scatter correction in computed tomography of 
industrial parts. 

[0002] Scatter is a deflection in radiation (suitable for imaging, such as X-rays) 
direction caused by certain interactions of the radiation within a target material, also 
referred to as an object. The phenomenon is significant in computed tomography of 
industrial parts because the materials, often metallic in nature, scatter x-rays to a 
greater degree, and is especially deleterious for three-dimensional, volumetric 
computed tomography where the entire object is irradiated by a cone beam of x-rays. 
This spatially-varying background adds to the true signal and can produce pronounced 
artifacts when the three-dimensional image of the object is mathematically 
reconstructed. 

[0003] The primary measurement data in computed tomography are sets of x-ray 
projections taken from various angles with respect to the object. In what follows, it is 
assumed that a complete computed tomography system having a computer is 
available, that the distribution of x-ray intensity in the various projection views has 
been detected, measured, and digitized in some manner known to the art to obtain raw 
projection data, which numerically represent the projection view of the object. These 
arrays of numbers representing the various projection views are accessible for 
numerical operation by the computer. It is further assumed that the projection arrays 
are afterwards combined and processed according to the known methods of computed 
tomography (CT) to produce a two-dimensional (2D) or a three-dimensional (3D) x- 
ray attenuation map, or a 2D or 3D representation, of the object. Typically generation 
of a single 2D image in a single reconstruction step is referred to as planar CT and 
generation of a 3D image in a single reconstruction step is referred to as volumetric or 
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cone-beam CT. Further a series of 2D images with appropriate step between adjacent 
slide locations is comparable to the 3D image set produced in a single cone-beam 
reconstruction. Typically, planar CT systems utilize a linear detector array and cone- 
beam CT systems utilize an area detector array. 

[0004] Hereinafter, reference to a 3D image will include by implication reference 
to a 2D image as a subset of a 3D image. Further, the computer is configured to 
provide the means for the reconstruction of and analysis of a voxellized representation 
of the object. 

[0005] It is known that the image artifacts caused by scattered x-rays falling on the 
various projection views can be corrected if the fraction of total signal at each point of 
every projection caused by scatter is estimated and then digitally subtracted before the 
projections are combined in the image reconstruction step. 

[0006] Until now, approaches for estimating this scattered component include 
making ancillary measurements using a series of x-ray blocking slits of varying width 
placed between the object and the x-ray detector. The rationale is that the scattered 
signal, being incident from a range of directions, can be estimated by extrapolating the 
series of slit measurements to zero width. However, such a method requires extensive 
added hardware and provides only a coarse grid of scatter estimates. More 
importantly, this approach has proven experimentally difficult and unable to provide 
accurate scatter estimates. 

[0007] A different approach has involved calculation of the scattered signal from 
physical first principles using prior knowledge of the object geometry. Accurate 
scatter estimates may be possible in this way using Monte Carlo radiation transport 
computer codes. However such estimates are calculation heavy and consume a large 
amount of processor time, a requirement which is prohibitive where a variety of 
different complex shape are to be imaged, as is the case in industrial imaging. 

[0008] It would therefore be desirable to have methods and systems that provide 
substantially accurate scatter correction estimates, and provide advantage in terms of 
computation time. 
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BRIEF DESCRIPTION 

[0009] Briefly, in accordance with one embodiment of the present invention, a 
method of correcting scatter includes obtaining a voxellized representation of a 3D 
image of an object from a plurality of projection data. A single scatter profile for the 
object is calculated using the voxellized representation of the 3D image of the object. 
A total scatter profile for the object is determined using the single scatter profile and 
an adjustment factor and the projection data is corrected using the total scatter profile 
to obtain a scatter corrected projection data. 

[0010] According to another embodiment, a beam hardening correction method 
includes simulating a number of attenuation data for an x-ray spectrum, at least one 
object material, and a detector spectral response. A function is fitted to the attenuation 
data to obtain an attenuation curve. A number of projection data for an object are 
corrected using the attenuation curve to obtain a number of beam hardening corrected 
projection data. A corrected image of the object is reconstructed from the beam 
hardening corrected image data. 

[0011] According to another embodiment, an imaging system for correcting scatter 
in an image of an object includes at least one radiation source adapted to expose the 
object to x-rays. A detector arrangement is disposed with respect to radiation source 
to receive x-rays passing from said radiation source through the object and a computer 
system is coupled to the detector arrangement. The computer system is further 
configured to acquire a number of projection data from the detector arrangement and 
generate a 3D image from the projection data. The computer system further generates 
a voxellized representation of the 3D image of the object, calculates a single scatter 
profile for the object using the voxellized representation of the 3D image of the 
object, determines the total scatter profile for the object using the single scatter profile 
and an adjustment factor and corrects the projection data using the total scatter profile 
to obtain, respectively, at least one of a scatter corrected projection data and a scatter 
corrected 3D image. 
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DRAWINGS 

[0012] These and other features, aspects, and advantages of the present invention 
will become better understood when the following detailed description is read with 
reference to the accompanying drawings in which like characters represent like parts 
throughout the drawings, wherein: 

FIG. 1 is a schematic illustration of an imaging system according to an 
embodiment of the invention; 

FIG. 2 illustrates a method for scatter correction according to an 
embodiment of the invention; 

FIG. 3 illustrates combination of voxels to form larger voxels; 

FIG. 4 illustrates removal of voxels having a CT number lower than a 
threshold value; 

FIG. 5 shows an exemplary x-ray attenuation curve; and 

FIG. 6 illustrates an exemplary beam hardening curve. 

DETAILED DESCRIPTION 

[0013] FIG. 1 shows an imaging system 50 for correcting scatter in an image of an 
object 20, according to an embodiment of the invention. At least one radiation source 
40 is adapted to expose the object 20 to radiation suitable for imaging, for example, x- 
rays. A detector arrangement 42 is disposed with respect to the radiation source 40 to 
receive x-rays passing from the radiation source 40 through the object 20. The x-rays 
passing through the object project a two-dimensional impression, referred to as raw 
projection data 10, on the detector arrangement 42. A computer system 44 is coupled 
to the detector arrangement 42 and configured to acquire a set of raw projection data 
10 from the detector arrangement 42, to generate a three dimensional (3D) image 12 
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from the raw projection data 10. The projection data 10 is comprised of a number of 
projections from different views, depending upon the application, and in some cases 
this number varies between 180 - 2500. However the number of projection views 
comprised in the projection data 10 is not restrictive on the present embodiment. The 
functionality of the computer system 44 is illustrated, for example, in FIG. 2. One 
skilled in the art will appreciate here that the raw projection data 10 may pertain to 
planar raw data or to cone beam raw data (the imaging system being a planar CT 
system or a volumetric CT system, respectively), and in either case, the computer 
system 40 is suitably configured to generate a 3D image 12 of the object 20. The 
computer system 44 is further configured to generate a voxellized representation 14 of 
the 3D image 12 of the object 20 and calculate a single scatter profile 13 for the object 
20 using the voxellized representation 14 of the object 20, as indicated in FIG. 2, for 
example. A total scatter profile 15 for the object 20 is obtained using the single 
scatter profile 13 and an adjustment factor, which relates the single scatter profile to a 
multiple scatter profile for the object 20. The computer system 44 is further 
configured to calculate the total scatter profile as above, and to obtain a scatter 
corrected 3D image 16. 

[0014] As used herein, "adapted to", "configured" and the like refer to mechanical 
or structural connections between elements to allow the elements to cooperate to 
provide a described effect; these terms also refer to operation capabilities of electrical 
elements such as analog or digital computers or application specific devices such as an 
application specific integrated circuit ASIC that are programmed to perform a sequel 
to provide an output in response to given input signals. 

[0015] According to a more specific embodiment of the disclosure, as illustrated 
by FIG. 2, a method of correcting scatter in industrial CT includes obtaining the 
voxellized representation 14 of the 3D image 12 of the object 20. The 3D image 12 is 
reconstructed from the raw projection data 10. A single scatter profile 13 for the 
object is calculated from an analytical method 1 1 using the voxellized representation 
14 of the object 20 as follows: primary beam intensity can be calculated for each 
voxel 22, 24 by calculating the attenuation between the radiation source 40 and the 
voxel 22, 24. (The voxels are indicated in FIG. 3, for example.) Total scattering x- 
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ray intensity for each voxel can be calculated using the scatter cross-section and 
material composition and density for the respective voxel. The single scatter profile 
for each detector pixel 38 is then calculated analytically using the scatter angular 
distribution, detector pixel solid angle, scattering x-ray attenuation and the detector 
response functions. The single scatter profile, so obtained, is used to determine a total 
scatter profile 15, which can be treated as a function of the single scatter profile. 

[0016] The total scatter profile is the sum of the single scatter and multiple scatter 
profiles. Various techniques, for example Monte Carlo, can be utilized to demonstrate 
that the multiple scatter profile is similar in shape to the single scatter profile for most 
industrial cases. Because of this similarity, the multiple scatter profile can be 
approximated as a product of the single scatter profile and an adjustment factor. The 
adjustment factor may be a constant or a variable function, and is arrived at by, among 
other techniques, three different approaches: (1) optimizing image quality of the 
object 20, the primary step proposed herein, (2) by calibration using outside 
collimated detectors, or (3) by deriving a total attenuation estimate for the entire beam 
traversing the entire object. In calibration processing (approach (2)), some outside 
detectors are blocked from receiving any primary beam from the x-ray source 40. The 
recorded signals from these detectors are due to the scatter contribution only. By 
comparing the difference between the actual measured scatter intensity and the 
calculated single scatter contribution for these detectors, the adjustable factor can be 
easily determined. In approach (3), the total attenuation of all rays is accumulated to 
give a single global approximate scatter factor, including attenuation suffered by 
scattered X rays as they exit the object. Optimizing image quality (or approach (1)) is 
flexible and accurate and typically requires comparatively less effort and time than the 
other approaches. In this approach the experimentally obtained projection data 10 is 
corrected using the total scatter profile, and a scatter corrected 3D image 16 of the 
object 20 is obtained. The correction is driven by attempting to obtain a relatively 
uniform background, and the degree of uniformity that is sought for obtaining the 
relatively uniform background is a selectable parameter depending upon a specific 
application. 
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[0017] According to a more specific embodiment, the total scatter is subtracted 
from the projection data 10 to obtain a scatter corrected projection data 18, which then 
is used to reconstruct 3D image 16. This is achieved as follows: the total scatter is 
represented as a sum of single scatter and multiple scatter, and the multiple scatter is 
represented as a product of the single scatter and an adjustment factor. The total 
scatter profile is subtracted from projection data 10, to obtain a scatter corrected 3D 
image 16. 

[0018] According to a more particular embodiment of approach 1, the adjustment 
factor is determined iteratively. For example, the corrected 3D image 16 of the object 
20 is analyzed. For example, optimizing image quality is used to analyze the 
corrected 3D image 16. For every variation of the adjustment factor, the corrected 3D 
image 16 is checked to determine to the degree expected improvements are obtained. 
Improvements might include establishing uniformity of cross sectional opacity in' 
cases where such is to be expected and/or reducing artifacts which are 
characteristically evident in CT images due to such effects, including rounding of 
edges and streaks aligned object features such as extended straight edges. If a clear 
improvement is not obtained, the adjustment factor is refined, and the total scatter 
profile is recalculated using the single scatter profile and the adjustment factor. The 
projection data 10 is corrected using a refined total scatter profile to obtain the scatter 
corrected 3D image 16. These steps are repeated until a satisfactory corrected image 
16 is obtained. 

[0019] The scatter correction method, as discussed, advantageously allows for 
incorporating other efficiency increasing techniques. According to an embodiment 
illustrated by FIG. 3, a voxellized representation 14 of the object 20 is obtained by 
combining a number of voxels 22 to form respective large voxels 24. In a related 
embodiment, a number of voxels 22 may be reconfigured to form the large voxels 24. 
This combination or reconfiguration reduces the number of voxels and results in 
speeding up calculations. This step is expected to only minimally impact the accuracy 
of the resultant calculated scatter profiles due to the fact that scattering is, by its 
nature, a low frequency process which is dependent more upon the global rather than 
the local distribution of material. According to another embodiment, an efficiency 
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increasing technique comprises correcting selected projection views from the 
projection data 10, and then interpolating the scatter correction of these projection 
views, to estimate the scatter corrections for the intervening views. For example, 
scatter profiles for thirty six equally spaced projections at ten degree intervals are 
generated, thereby covering projection views spread over 360 degrees. Thereafter, 
scatter profile for the remaining views is interpolated from the scatter profile of the 
thirty six projection views, which advantageously speeds up the scatter correction 
process. 

[0020] According to another embodiment illustrated in FIG. 4, obtaining the 
voxellized representation 14 of the object 20 includes determining a threshold for the 
3D image 12. A CT number, which is a measure for the material being examined, is 
determined for each of the voxels 22, 26, 28, and compared to the threshold value. 
Based on the comparison, some of the voxels 26, 28 are removed to simplify the 
voxellized representation of the object. For example, the threshold value may be set 
at the CT number for air, and accordingly, all the voxels that contain only air may be 
removed, because these voxels do not significantly contribute to any attenuation or 
scatter in the x-rays. Thus, the voxels having an insignificant contribution to the 
object image and the scatter are eliminated, reducing the calculation time. 

[0021] It will be appreciated that the computer system 44 is appropriately 
configured to execute the methods and techniques for scatter correction, as discussed. 
Briefly, the imaging system 50 includes a computer system 44 adapted to receive a 
number of raw projection data 10 from the detector arrangement 42, and further 
adapted to correct the projection data 10 using the total scatter profile to generate a 
corrected 3D image 16, as discussed. The computer system 44 is configured for 
executing various numerical operations. The numerical operations include calculation 
the single scatter profile, the total scatter profile for the object 20 by multiplying the 
single scatter profile by the adjustment factor, the corrected projection data 10 
obtained by subtracting the total scatter profile from the projection data, and 
reconstructing the scatter corrected 3D image 16. The numerical operations also 
include obtaining the voxellized representation 14 of the object 20 by combining a 
number of voxels 22 to form large voxels 24, removing voxels 26, 28 based on 
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comparison of CT numbers to a threshold value, to simplify the voxellized 
representation 14 of the object. The computer system 44 is further configured to 
iteratively calculate the corrected image 16 of the object by analyzing the corrected 
image 16, refining the adjustment factor based on the analysis, recalculating the total 
scatter profile for the object 20 using the single scatter profile and the adjustment 
factor and correcting projection data 10 using the total scatter profile to provide a 
corrected projection data 18 set which is then used to reconstruct the scatter corrected 
3D image data 16. 

[0022] In addition to scatter, CT images are degraded by "beam hardening" effects 
that result from the use of broadband or poly-energetic x-ray sources. As shown for 
example in FIG. 5, by the absorption coefficient versus incident energy graph 46, 
lower energy x-ray photons are preferentially absorbed. Consequently, longer path 
lengths through an object will exhibit lower attenuation. Exemplary beam hardening 
curves depicting ideal intensity attenuation versus path length 47 and actual intensity 
attenuation versus path length 48 are shown in FIG. 6. The beam hardening curve 
depends on the x-ray spectrum, the object material, and the detector spectral response. 
According to another embodiment of the invention, a beam hardening correction 
method is provided. The method includes simulating a set of attenuation data for an 
x-ray spectrum, at least one object material and a detector spectral response. A 
function is fitted to the attenuation data to obtain an attenuation curve. For example, a 
polynomial fit, such as a third order polynomial, can be employed. A set of projection 
data 10 for an object 20 are corrected using the attenuation curve to obtain a set of 
beam hardening corrected projection data. A corrected image of the object is 
reconstructed from the beam hardening corrected projection data. 

[0023] Though the scatter and beam hardening correction methods and systems 
described herein have been explained with respect to scatter as encountered in planar 
(2D) or cone-beam (3D) CT systems, one skilled in the art will appreciate the method 
applies to all radiation imaging techniques, for example, radiographic imaging, with 
linear arrays or area detectors, among others. 
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[0024] Although only certain features of the invention have been illustrated and 
described herein, many modifications and changes will occur to those skilled in the 
art. It is, therefore, to be understood that the appended claims are intended to cover all 
such modifications and changes as fall within the true spirit of the invention. 
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